Spin-orbit-enhanced Wigner localization in quantum dots 
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We investigate quantum dots with Rashba spin-orbit coupling in the strongly-correlated regime. 
We show that the presence of the Rashba interaction enhances the Wigner localization in these 
systems, making it achievable for higher densities than those at which it is observed in Rashba-free 
quantum dots. Recurring shapes in the pair distribution functions of the yrast spectrum, which 
might be associated with rotational and vibrational modes, are also reported. 
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I. INTRODUCTION 

Electron localization in finite-size nanostructures, rem- 
iniscent of the Wigner crystallization of the bulk electron 
gasi, occurs when the density of the system becomes 
low enough so that the Coulomb repulsion dominates 
over the motion associated with the kinetic energy of 
the electrons. As a consequence, the electrons "crystal- 
lize", acting as classical charges seeking for the equilib- 
rium positions that minimize their total energy. From 
the theoretical point of view, this phenomenon has been 
extensively investigated in e.g. quantum dots (QDs) 2 ~ 5 
(see also the review Ref. 6), quantum wiresi, and quan- 
tum rings^, and analogue Wigner-localized states have 
also been observed in other quantum systems such as 
e.g. ultracold atomic gases, ultracold trapped ions, com- 
plex plasmas, or cuprate chain compounds 9 . On the ex- 
perimental side, signatures of electron localization have 
been found in e.g. quantum dotsi^, quantum wires^i, or 
carbon nanotubesi^. The study of this phenomenon is 
not only interesting from the purely fundamental point 
of view, but it has also been shown that Wigner-localized 
systems can be useful to design e.g. chains of spatially 
separated quantum bitsi^ or quantum hard drives^. 

In experiments with semiconductor quantum dots and 
wires, however, difficulties are encountered when lower- 
ing the density down to the regimes at which Wigner 
localization is achieved due to disorder effect o 10 ' 11 . The 
search for more clear and direct experimental evidence of 
Wigner localization in such systems, as well as for ways 
that allow its observation at more easily accessible den- 
sities, is therefore an interesting matter of research. 

In the context of semiconductor nanostructures, the 
study of the effects of the so-called Dresselhaus 14 (DSOI) 
and Rashba^ spin-orbit interactions (RSOI) has also 
attracted great interest in the last years. The origin 
of these interactions is due to the existence of internal 
asymmetries that give rise to net electric fields within 
the heterostructures. In the reference frame of the elec- 
trons, these electric fields transform into effective mag- 
netic fields, which couple to the electronic spins. In the 
case of the Rashba spin-orbit coupling, which is due to 
the asymmetry of the potential well confining the con- 
duction electrons within the heterostructure, special at- 



tention has been paid since it was shown that its strength 
can be experimentally tuned within large ranges of values 
by means of the application of external electric fields^ or 
by changing the electron density— , making of this mech- 
anism a very promising tool for potential applications in 
e.g. spintronics or quantum computation since it pro- 
vides a natural way to control and manipulate the elec- 
tronic spir>i£. 

Important effects due to the spin-orbit (SO) interac- 
tion in semiconductor nanostructures have been reported 
on, e.g., the addition energies, the electron g- factor, 
the magnetoconductivity, the spin textures, or the spin- 
relaxation properties^. Also, spin-orbit-induced electron 
localization has been observed in Coulomb-free closed 
loops of quantum wires with Rashba SO coupling, as 
well as in quantum rings with both RSOI and DSOI 20 . 
It must be pointed out, however, that this localization is 
not of the Wigner-type since the electrons are considered 
to be non-interacting, but is entirely due to the analo- 
gies and the interplay of the spin-orbit coupling with 
external magnetic fields, which gives rise to Aharonov- 
Bohm cages and periodic trapping potentials that cause 
the localization of the particles^ -. As a matter of fact, 
the electron-electron interaction has been taken into ac- 
count in a rather small fraction of the many studies on 
the spin-orbit coupling in semiconductor nanostructures, 
and only a few have treated it without approximations 
or have considered its interplay with the SO terms 2 ^. 

In this Letter, we investigate strongly-correlated InAs 
quantum dots with RSOI and we show that Wigner local- 
ization is largely enhanced by the presence of the Rashba 
coupling, making this phenomenon potentially observ- 
able at sensibly higher electronic densities than those 
at which it is achieved in spin-orbit-free QDs. Since 
the Dresselhaus contribution is expected to be small in 
narrow-bandgap materials such as InAs 2 ^, we neglect it 
in the present study. 



II. MODEL AND METHOD 

We consider a two-dimensional InAs quantum dot with 
parabolic confinement in the xy— plane given by V(r) = 



(1/2V 2 



y ). The Rashba spin-orbit coupling is de- 
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scribed by the usual terrain 

Hr = [p y a x - p x <T y ] , (1) 

m*m e 

where m*m e is the electron effective mass, with m* = 
0.023, pi is its linear momentum in the i— direction, ui is 
the corresponding Pauli matrix, and kR determines the 
strength of the interaction. 

The QD is considered to have N electrons interact- 
ing with each other through the Coulomb interaction. 
Expressing the energies and lengths in units of Tvuj and 

l u = \J m J^ n J respectively, the full Hamiltonian of the 

system, in second quantization form, reads 

a/3 

+ k R [p y a x - p x a y ] | (3^ a^ap 

+ \ ^ $ I I r - r' I ' 7 ' ^ a tx a l a "i a & > ( 2 ) 

a/3f8 ' ' 

where aL (p = a, /3, 7, S) are the usual creation 
and annihilation operators for the single-particle state 
and fcfl = l^kjf is a dimensionless parameter de- 
termining the effective spin-orbit interaction strength. 
The prefactor in the Coulomb term is given by g — 
e 2 / [A-neQerluhuS) = l^/a* B , where a* B is the effective Bohr 
radius, which for InAs is approximately 34 nm. The 
confinement frequency is related to the average electron 
density n through the relation to 2 = fi 2 /(m 2 a* B VNr 3 ), 3 
where r s = 1/ ^/im is the Wigner-Seitz radius. 

Using the configuration-interaction method, we nu- 
merically find the many-particle eigenstates of the sys- 
tem, expanded in a basis of Slater determinants made up 
of 2D harmonic-oscillator single-particle orbitals \fi) = 
\n, m, s), where n, m and s =t,| are, respectively, the 
radial, azimuthal, and spin quantum numbers. Defining 
the equivalent set n± = (2n + |m| ± m)/2, one has for 
the matrix elements of the Rashba term: 

(n', m',t| Hr I n, m,f) = (n', m',i\ H R \ n, m,I) = 0, 
(n+',nJ, t| H R I n +) n_4) = ~k R x 

(n+',nJ,i\ H R I n+,n_,t) = k R x 

](3) 

The (in principle) infinitely large basis space is trun- 
cated, taken to be large enough to ensure converged so- 
lutions, which in our calculations are achieved with 15 
harmonic oscillator shells. As the total angular momen- 
tum J z = L z + S z represents a good quantum number, 
the diagonalization is in practice done separately for each 
considered J z . 

A standard value for the confinement energy is of the 
order of Tiuj ~ 5 meV, with associated effective Coulomb 



strength g = 0.75 = go. We have studied the interval 
9o < 9 < 3(?o , corresponding to 5 meV < fko < 0.56 meV. 
Finally, regarding the RSOI term, we have considered 
the range 0.263 < kR < 2.63, and we have found that 
the spin-orbit-induced localization effects become clearly 
observable for strengths of the Rashba coupling constant 
h 2 kR/m > 8.5 x 10~ u cV-m. While being slightly more 
than two times larger than the ones experimentally ac- 
cessible nowaday o 17 ! 23 , such values should be reachable 
in the near future considering the ongoing research on 
new materials with larger SO coupling 

III. RESULTS 

When the RSOI is neglected (kR = 0), we find that 
the electrons start to localize only for strengths of the 
Coulomb interaction larger than ~ 5go, corresponding to 
r s > AAa B , as expected 3 -^. For weaker interactions, we 
only observe a slight deformation of the electron cloud 
as a whole. Also, when the localization occurs, there is a 
depletion of the electron density in the center of the trap 
due to the repulsion between the particles. 

This scenario, however, sensibly changes when the 
Rashba term is taken into account. In Figs. 1 and 2 
we show the pair distribution function 

PDF{y,y') = $\t)${t')${t')${t)) , (4) 

which gives the probability of finding a particle at the 
position r provided that another one is at r', for the 3- 
and 4-electron quantum dots for different strengths of the 
effective spin-orbit and Coulomb interactions. One can 
clearly see that the inclusion of the Rashba term, and 
especially for large values of kR, significantly enhances 
the electron localization. Indeed, in both cases it is al- 
ready visible for kR = 1.84 and g = 3gch corresponding 
to Tiuj = 0.56 meV and h 2 k R /m = 8.5 x 10 -11 eV-m, with 
the electrons distributing themselves on the vertices of an 
equilateral triangle and of a square, respectively, analo- 
gously to what is observed in quantum dots submitted 
to perpendicular magnetic fields 2 ' 10 . This corresponds 
to Wigner-Seitz radii of r s = 2.5a* B for N = 3 and to 
r s = 2Aa* B for N — 4, i.e. almost half the values at 
which localization is observed in spin-orbit-free quantum 
dots. A further increase of the RSOI while keeping the 
confinement strength fixed leads to a more pronounced 
localization of the electrons, as can be seen from the fig- 
ures when kR = 2.63. 

We have also investigated a dot with six electrons, al- 
though in this case an optimal convergence of the numer- 
ical simulations has not been possible to achieve. Nev- 
ertheless, we have observed that the behavior of the sys- 
tem is qualitatively similar to the cases with TV = 3 and 
N = 4, with localized configurations appearing at simi- 
lar effective Coulomb strengths with r s ~ 2.2o^, made 
up of five electrons forming a pentagon with the sixth 
one placed at the center of the dot. Finally, we also 
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want to mention that in some papers electron localiza- 
tion has been addressed making use of quantitative suit- 
able criteria^S. In our study, however, we have limited 
ourselves to the rather qualitative discussion presented 
above. 
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FIG. 1: Pair distribution functions for the 3-electron quantum 
dot and for different strengths of the effective Rashba (increas- 
ing from top to bottom) and Coulomb (increasing from left 
to right) interactions. The star symbol indicates the position 
of the reference electron. 



The enhancement of the Wigner localization due to the 
Rashba spin-orbit coupling can be qualitatively explained 
already from the non-interacting-electron picture. In- 
deed, we have observed that the spatial extension of the 
single-particle density profiles is reduced when the RSOI 
term is taken into account as compared to the harmonic 
oscillator case with kji = 0. This is in close analogy 
with the Fock-Darwin single-particle states of quantum 
dots under magnetic fields 2 - and may be understood from 
the nature of the spin-orbit interaction described above. 
Therefore, it can be concluded that one of the effects 
of the Rashba coupling is to minimize the overlapping 
between the different electrons and thus to favour their 
spatial separation with respect to each other. 

It is also interesting to study the so-called yrast 
spectrumSi, i.e., the lowest energy level for a given to- 
tal angular momentum, which we plot in Fig. 3 for the 
three-electron dot in different parameter regimes. This 
again shows very close analogy with what is observed in 
quantum dots under perpendicular magnetic fields^, the 
main difference being in the low J z range. In the Rashba 
case, low \j z \ single-particle states are almost degenerate 
and therefore the spectrum is rather flat for low values 
of J z . On the contrary, under an applied magnetic field, 



FIG. 2: Same as Fig. 1 for the quantum dot with N = 4. 



a preferred spin direction is selected and low J z states, 
realized using s.p. states of opposite angular momentum, 
are energetically unfavourable. This explains the overall 
parabolic shape of the yrast profiles shown in^, different 
from the almost degenerate low J z levels in the Rashba 
case. For higher angular momenta, however, the mag- 
netic and Rashba yrast spectra become more similar: in 
both cases the yrast lines show oscillations with period 
N, which become more pronounced as the Coulomb re- 
pulsion increases. The presence of these oscillations is 
explained in Ref. 2 by the vanishing of the exchange 
energy term for configuration with electrons occupying 
adjacent j z single-particle levels, which occur only for ev- 
ery third unit of angular momentum. Moreover, for the 
strong magnetic field case -and thus for high J z - it can 
be seen that the electrons behave as a rigid, molecule-like 
configuration 2 ^. A semiclassical interpretation of the 
yrast spectrum is then possible, with the local minima 
corresponding to purely rotational modes and with the 
intermediate states containing at least one vibrational 
quantum 29 . 

Finally, a qualitative signature of the above-mentioned 
molecular behaviour in the Rashba-interacting case can 
be observed from Fig. 4, where we show the pair distri- 
bution functions corresponding to the Wigner-localized 
case for different values of the total angular momen- 
tum. One can recognize recurring shapes (J z — 7.5ft 
and J z = 10.5ft, J z = 8.5ft and J z = 11.5ft, J z = 9.5ft 
and J z =12.5), which might be associated with specific 
rotational and vibrational states for localized electrons 
in quantum dots. Indeed, the standard triangular shapes 
are found at the local minima, consistently with the hy- 
pothesis that they correspond to rotational modes, where 
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J z (h) 

FIG. 3: Yrast spectrum for the 3-electron quantum dot with 
g = 3go and kn — (blue diamonds), g = and kn = 2.63 
(green squares), g = go and kn = 2.63 (red triangles), and 
g = 3(?o and kn = 2.63 (black circles). The energies have 
been shifted to have the same value at J z = 0.57i for a better 
comparison. 
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FIG. 4: Pair distribution functions for different values of J z 
for the three-electron dot in the Wigner-localized regime (g = 
3(7o and kn = 2.63). 



no deformation is involved. Other shapes also appear pe- 



riodically, and they look consistent with the vibrations of 
a triangular molecule (stretching and shrinking the tri- 
angle basis, for example, similarly as patterns previously 
observed 28 ). We think that the appearance of these de- 
formed pair distribution functions matching the period- 
icity of the yrast spectrum is a strong hint of a "rigid- 
molecule" behavior of the electrons. 



IV. CONCLUSIONS 

In conclusion, we have found that the presence of 
Rashba spin-orbit coupling in strongly-correlated quan- 
tum dots largely enhances the Wigner localization in 
these systems and makes it observable at higher densi- 
ties than in the spin-orbit-free case. We report several 
analogies with quantum dots under magnetic fields, and 
in particular we show that the pair distribution func- 
tions corresponding to the states of the yrast spectrum 
of the Wigner-localized dot present recurring shapes that 
one can associate with rotational and vibrational modes. 
Our results thus point out the importance of the exper- 
imental investigation of Wigner localization in materials 
with large spin-orbit coupling. Although not modeled 
here, any disorder in the material can be expected to 
contribute to further enhance the localization 10 ' 11 . The 
effects of the Dresselhaus spin-orbit interaction on the 
electron localization, and in particular its interplay with 
the Rashba term, is also an interesting subject that will 
be investigated in future works. 
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